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Abstract. With the acceleration of lead nuclei in the LHC, heavy- ion physics will 
enter a new energy domain. One of the main novelties introduced by the 30-fold 
£h \ energy-jump from RHIC to the LHC is the abundant heavy-quark production. After 

1 ^ ' discussing a few examples of physics issues that can be addressed using heavy quarks, 

we present a selection of results on the expected experimental capability of ALICE, 
| the dedicated heavy-ion experiment at the LHC, in the open-heavy-flavour sector. 
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1. Introduction: novel aspects of heavy-ion physics at the LHC 

in ■ 

The nucleon-nucleon c.m.s. energy for Pb-Pb collisions at the LHC, y/s NN = 5.5 TeV, 
will exceed that available at RHIC by a factor about 30, opening up a new domain for 
the study of strongly-interacting matter in conditions of high temperature and energy 
density (QCD medium). Particle production at the LHC will present quantitative and 
^ 1 qualitative new features, as discussed in the following. 

Hard processes should contribute significantly to the total cross section. The 
^ ■ mechanism of energy loss due to medium-induced gluon radiation allows to use the 
energetic partons produced in initial hard-scattering processes as probes to collect 
information on the opacity and density of the medium itself. At the LHC, the set 
of available probes will be extended both quantitatively and qualitatively. In fact, hard 
(light)quarks and gluons will be produced with high rates up to very large transverse 
momentum (p t ). Additionally, charm and beauty quarks, which, due to their masses, 
would show different attenuation patterns (see section |2J), will become available for 
detailed measurements, since their production cross sections are expected to increase by 
factors 10 and 100, respectively, from RHIC to the LHC [1]. 

High-density parton distributions are expected to dominate particle production. 
The density of low momentum-fraction, x, gluons in the two colliding nuclei is expected 
to be close to saturation of the available phase space, so as to produce significant 
recombination effects. As an example, at central rapidity, low-p t c (b) quarks will 
be produced by partons with x\ ~ x 2 Z 2m c (mb)/ v / % ^ 2.4(10) GeV/5500 GeV ~ 
4(16) x 10 -4 . In section 12 we will discuss how gluon fusions in this x region are expected 
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Table 1. Expected QQ yields at the LHC, from NLO pQCD (parameters in the 
text) [5]. For p-Pb and Pb-Pb, PDF shadowing is included and N co \\ scaling is applied. 



colliding system a/s N n centrality iV cc /event iV bb /event 



pp 14 TeV minimum bias 0.16 0.0072 

p-Pb 8.8 TeV minimum bias 0.78 0.029 

Pb-Pb 5.5 TeV central (0-5% a tot ) 115 4.6 



to affect cc (and, to a lower extent, bb) production, not only in p-Pb and Pb-Pb 
collisions (nuclear shadowing), but possibly even in pp collisions. 

2. Heavy-quark phenomenology from pp to nucleus— nucleus collisions 

Heavy-quark pairs, QQ, are produced in partonic scatterings with large virtuality 
(momentum transfer) Q > 2mq. Therefore, the production cross sections in 
nucleon-nucleon (NN) collisions can be calculated in the framework of collinearly 
factorized perturbative QCD (pQCD). The differential QQ cross section is written as: 



which depends on the renormalization scale currently, calculations are performed 
up to next-to-leading order (NLO), 0(a^). The nucleon Parton Distribution Function 
(PDF) for the parton of type i at momentum fraction X\ and factorization scale fi-p, 
which can be interpreted as the virtuality of the hard process, is denoted by /i(xi,/i|). 

The expected yields in pp collisions at y/s = 14 TeV are reported in the first line 
of Table ^ These numbers are obtained at NLO using the MNR program [2] with 
m c = 1.2 GeV and \ip = = 2 m c for charm and m b = 4.75 GeV and fip = fiR = 
for beauty; the PDF set is CTEQ 4M [3]. The predicted yields have large uncertainties, 
of about a factor 2, estimated by varying the values of the masses and of the scales 
(much smaller uncertainties, ~ 20%, arise from the indetermination in the PDFs) [4,5]. 
The theoretical uncertainty band for the D-meson cross section as a function of p t will 
be shown in sectional along with the expected sensitivity of the ALICE experiment [6]. 

As aforementioned, in the x range relevant for QQ production, the PDFs will be 
close to phase-space saturation and, already in the case of pp collisions, there will be 
important gluon-fusion effects (gg — > g). These can be accounted for in the PDF scale- 
evolution equations by adding to the standard linear DGLAP term a negative nonlinear 
(quadratic) term (see [7] and references therein): 





df g (x, Q 2 )/d\ogQ 2 = [DGLAP term of 0(/ g )] - [term of C(/ g 2 )] . (2) 
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Figure 1. Left: enhancement due to nonlinear gluon evolution for c quarks and D 
mesons in pp collisions at y/s — 14 TeV [9]. Centre: modification of the gluon PDF in 
a Pb nucleus at Q 2 ~ 4m;?. Right: corresponding R^ A in p-Pb at ^/s NN = 8.8 TeV. 



The nonlinear term, currently calculated only at LO, 'slows down' the Q 2 evolution at 
given x. It has been shown [7] that, for x < 10~ 2 , it allows to have a higher gluon density 
at small Q 2 ( ^ 10 GeV 2 ), with respect to that obtained with DGLAP terms only, and to 
maintain at the same time a good fit of the proton structure function data from HERA. 
A higher gluon PDF would imply an enhancement, w.r.t. to DGLAP-based calculations, 
of cc production at low p t at LHC energy [8]. Figured] (left) shows, as a function of 
Pt, the enhancement at the c-quark and at the D- meson level, for m c = 1.2 GeV and 
Mf = Mr = Q 2 = 4m t 2 c = 4 {m 2 + p 2 ) [9]. Here, as well as for the other results 
presented in the following, the hadronization of heavy quarks is performed using the 
string fragmentation model implemented in PYTHIA [10]. The enhancement survives 
fragmentation and it is of about 30% for D-meson p t — > 0, even in this 'pessimistic' case 
where relatively-large Q 2 values are considered (it should be noted, however, that this 
is a LO result and the effect might be smaller at NLO). In section El we will discuss a 
possible strategy to detect the enhancement in pp collisions at the LHC. 

For hard processes, in the absence of nuclear and medium effects, a 
nucleus-nucleus (or p-nucleus) collision would behave as a superposition of independent 
NN collisions. The charm and beauty differential yields would then scale from pp to AA 
(or pA) proportionally to the number iV C oii of inelastic NN collisions (binary scaling): 

d 2 ^AA( P A)/dptd2/ = N coll x d^/dptdy . (3) 

Binary scaling is, indeed, expected to break down due to both initial-state effects, such 
as nuclear shadowing of the PDFs, and final-state effects, such as parton energy loss in 
the medium formed in nucleus-nucleus collisions. As a consequence, as we will detail in 
the following, heavy quarks are important tools to probe and investigate these effects. 

In Table Q we report the cc and bb yields in p-Pb and Pb-Pb collisions calculated 
including in the NLO pQCD calculation the EKS98 parameterization [11] of the PDFs 
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Figure 2. Nuclear modification 
factors for charged hadrons [15], D 
mesons (m c = 1.2 GeV) and B mesons 
(m b = 4.8 GeV) [16] in central (0- 
10% cr tot ) Pb-Pb collisions relative 
to binary scaling from pp collisions, 
at y/s NN = 5.5 TeV. The D- meson 
result does not include feed-down 
from B — * D + X decays, which is, 
however, expected to be rather small 
(~ 5%) [18]. 



nuclear modification ff b (x,Q 2 )/ ff(x,Q 2 ), shown in Fig. [T] (centre) for Q 2 = 5 GeV 2 , 
and applying binary scaling [5]. The charm (beauty) cross-section reduction induced 
by shadowing is about 35% (20%) in Pb-Pb and 15% (10%) in p-Pb. There is a 
significant uncertainty on the strength of shadowing in the small-x region and some 
models predict much larger suppression than EKS98 (see [12] for a review). The 
comparison of QQ production in pp and p-Pb collisions (where final-state effects are not 
present) is regarded as a sensitive tool to probe nuclear PDFs at the LHC. The ratio of 
invariant-mass spectra of di-leptons from heavy-quark decays in p-Pb and pp collisions 
would measure the nuclear modification f^ h / f£ [12]. Another promising observable in 
this respect is the nuclear modification factor of the D-meson p t distribution: 

pA(AA) (ft) - ^ X dm » /dptdy ■ (4) 

We note that at the LHC, the pp, p-Pb and Pb-Pb runs will have different yj s NN 
values; however, pQCD can be used to extrapolate the measured cross sections between 
different energies [4] and, thus, calculate the -R p a(aa) ratios. In Fig. Q (right) we show the 
sensitivity of R^ A to different shadowing scenarios, obtained by varying the modification 
of the PDFs (shown for gluons in the central panel of the same figure). 

Experiments at RHIC have shown that the nuclear modification factor Raa is an 
effective tool for the study of the interaction of the hard partons with the medium 
produced in nucleus-nucleus collisions. Heavy-quark medium-induced quenching is one 
of the most captivating topics to be addressed in Pb-Pb collisions at the LHC, where 
both c and b quarks will be produced with high rates (see Table HJ). Due to the QCD 
nature of parton energy loss, quarks are predicted to lose less energy than gluons (that 
have a higher colour charge) and, in addition, the 'dead-cone effect' is expected to reduce 
the energy loss of massive quarks [13, 14]. Therefore, one should observe a pattern of 
gradually decreasing Raa suppression when going from gluon-originated light-flavour 
hadrons (h) to D and to B mesons: R AA £ i?° A < -Raa- m Fig- El we report recent 
estimates of these quantities for central Pb-Pb collisions at the LHC [15,16], obtained 
in the framework of a model [15] where energy loss is simulated in a parton-by-parton 
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approach combining the BDMPS 'quenching weights' [17] and a Glauber-model-based 
definition of the in-medium parton path length. For c and b quarks, the quenching 
weights were specifically calculated using the formalism developed in [14]. The BDMPS 
transport coefficient (a measure of the medium density) at the LHC was set to the value 
q = 100 GeV 2 /fm, estimated on the basis of the analysis of RHIC data performed in [15]. 
The results are plotted as bands that represent the theoretical uncertainty [15,16]. For 
charged hadrons, Raa is given for p t > 8 GeV because at lower p t the soft-particle 
component produced from the radiated gluons, which is not implemented in the model, 
might strongly contribute in shaping the nuclear modification factor. Whereas for D 
and B mesons Raa can be calculated down to pt — by assigning a 'thermal' transverse 
momentum (according to dN/dm t oc m t exp(— m t /T), T = 300 MeV) to the c and 
b quarks that lose most of their initial energy in the medium [16]. The presence of 
this thermalized component determines the rise of Raa at — ► 0. It should be 
mentioned here that, in the low p t region, the hadronization of heavy quarks in nucleus- 
nucleus collisions is likely to happen inside the medium via parton recombination [19], 
thus producing deviations from the simple pattern obtained from the aforementioned 
thermalization assumption. 

3. Looking for heavy quarks at the LHC: tools, techniques, performance 

Three experiments will participate in the LHC heavy-ion program: ALICE, the 
dedicated heavy-ion experiment [6]; CMS, with a strong heavy-ion program [20]; (most 
probably) ATLAS, which recently expressed interest in participating [21]. The three 
detectors have different features and design requirements, but all of them are expected to 
have excellent capabilities for heavy-flavour measurements. Their complementarity will 
provide a very broad coverage in terms of phase-space, decay channels and observables. 

Experimentally, the two key elements for a rich heavy-flavour program are: 
tracking/vertexing and particle identification (PID). 

Open charm and beauty mesons have typical life-times of few hundred microns 
(cr values are about 125-300 /mi for D mesons and 500 /im for B mesons) and the 
most direct detection strategy is the identification of single tracks or vertices that are 
displaced from the interaction vertex. The detector capability to perform this task is 
determined by the impact parameter§ (do) resolution. All experiments will be equipped 
with high position-resolution silicon-detector layers, including pixels, for precise tracking 
and impact parameter measurement also in the high-multiplicity environment of central 
Pb-Pb collisions. Tracking is done in the central (pseudo)rapidity region: \r]\ < 0.9 
for ALICE and \r]\ < 1.5 for CMS and ATLAS. In Fig. El we show the do resolution of 
ALICE and CMS||, along with the p t resolution, which is another important ingredient 
for heavy-flavour measurements (e.g. it determines the invariant mass resolution). The 

§ We define as impact parameter the distance of closest approach to the interaction vertex of the track 
projection in the plane transverse to the beam direction. 

|| ATLAS may perform similarly to ALICE, but systematic Pb-Pb studies are still in progress. 
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Figure 3. Transverse momentum (left) and transverse track impact parameter (right) 
resolutions for the ALICE and CMS detectors in Pb-Pb collisions. 



main difference between the two experiments is given by the magnetic field values: 
ALICE (0.4 T) has a very low p t cutoff of 0.2 GeV, while CMS (4 T) has a higher cutoff 
of 1 GeV but better p t resolution at high p t . The do resolutions are quite similar and 
better than 50 /im for p t > 1.5-3 GeV. 

Both lepton and hadron identification are important for heavy-flavour detection. D 
and B mesons have relatively-large branching ratios (BR) in the semi-leptonic channels, 
~ 10% to electrons and ~ 10% to muons, and inclusive cross-section measurements 
can be performed via single leptons or di-leptons. ALICE can identify electrons with 
p t > 1 GeV and < 0.9, via transition radiation and dE/dx measurements, and 
muons in the forward region, 2.5 < 77 < 4, which allows a very low p t cutoff of 
1 GeV. CMS and ATLAS have a broad pseudorapidity coverage for muons, |?7| < 2.4 
and \r]\ < 2.7, respectively, but they have a quite-high p t cutoff of ~ 4 GeV. Both 
CMS and ATLAS have high-resolution electro-magnetic calorimeters that can be used 
to identify electrons, although performance studies for heavy-ion collisions have not 
been carried out yet. Semi-leptonic inclusive measurements do not provide direct 
information on the D(B)-meson p t distribution, especially at low p t , because of the 
weak correlation between the lepton and meson momenta. Therefore, for charm in 
particular, the reconstruction of exclusive (hadronic) decays is preferable. In this case, 
in a high-multiplicity environment, hadron identification allows a more effective rejection 
of the large combinatorial background in the low-p t region. ALICE disposes of 7r/K/p 
separation via dE/dx and time-of-flight measurement for p < 3-4 GeV and \rj\ < 0.9. 

In the following we present results on the expected performance for the detection of 
D and B mesons in ALICE (CMS and ATLAS studies are described in [4,22]). The cc 
and bb cross sections from Table Q are used and the charged-particles rapidity density 
of a central (0-5% <7 tot ) Pb-Pb collision is assumed to be dN ch /dy ~ 6000. The results, 
based on realistic detector simulations, are given for samples of 10 9 pp events (about 9 
months of data taking) and of 10 7 central Pb-Pb events (about 1 month of data taking). 
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Charm reconstruction in ALICE. One of the most promising channels for open charm 
detection is the D° — > K~tt + decay (and charge conjugate) that has a BR of 3.8%. The 
expected yields (BRxdN/dy at y = 0), in central Pb-Pb (0-5% <x tot ) at ^/s^ = 5.5 TeV 
and in pp collisions at y/s = 14 TeV are 5.3 x 10 _1 and 7.5 x 10~ 4 per event, respectively. 

The main feature of this decay topology is the presence of two tracks with 
impact parameters d ~ 100 /mi. The detection strategy [23] to cope with the 
large combinatorial background from the underlying event is based on the selection 
of displaced-vertex topologies, i.e. two tracks with large impact parameters and good 
alignment between the D° momentum and flight-line, and on invariant-mass analysis 
to extract the signal yield. This strategy was optimized separately for pp and Pb-Pb 
collisions, as a function of the D° transverse momentum, and statistical and systematic 
errors were estimated [18]. 

Figure 0] (left) shows the expected sensitivity of ALICE for the measurement of 
the D° ^-differential cross section in pp collisions, compared to the pQCD calculations 
uncertainty that we mentioned in section El and in central (0-5% cr tot ) Pb-Pb collisions. 
The low-p t reach, provided by the moderate ALICE magnetic field (0.4 T) and the 
K/tt separation via time-of-flight, allows to address the issue of charm 'enhancement' 
in pp collisions due to nonlinear gluon evolution (section HJ). Here, the idea is that the 
effect — enhancement only at low p t — cannot be mimicked by standard DGLAP-based 
pQCD just by varying the input parameters. This is illustrated in Fig. 0] (from [9]), 
where the data-to-theory ratio is plotted versus the D° p t . The data are obtained with 
m c = 1.2 GeV and Q 2 = ^m 2 c and they include the enhancement (from Fig. (left)), 
while the theory results do not: only for m c < 1.1 GeV the theory can mimic the 
enhancement, but such small values are not supported by lower-y's measurements. 

The direct measurement of the D-meson p t distribution allows a good sensitivity to 
the suppression of the nuclear modification factor due to c-quark energy loss in Pb-Pb 
collisions. Figure (left) shows the estimated experimental errors on reported on 
a curve that lies in the middle of the D-meson band from Fig. |21 

Simulation studies to assess the performance for D° reconstruction in p-Pb are in 
progress and the results are expected to be similar to the pp case [18]. This might allow 
to measure nuclear shadowing via the R^ A ratio, as discussed in section |2l 

Beauty via single electrons in ALICE. The expected yield (BR x dN/dy at y = 0) for 
B -> e ± + X per central (0-5% a tot ) Pb-Pb collision at v ^ = 5.5 TeV is 9 x 10" 2 . 

The main sources of background electrons are: (a) decays of D mesons; (b) decays 
of light mesons (e.g. p and u); (c) conversions of photons in the beam pipe or in the 
inner detector layers and (d) pions misidentified as electrons. Given that electrons 
from beauty have average impact parameter do — 500 pm and a hard momentum 
spectrum, it is possible to obtain a high-purity sample with a strategy that relies on: 
electron identification with a combined dE/dx and transition radiation selection, which 
allows to reduce the pion contamination by a factor 10 4 ; impact parameter cut to reject 
misidentified pions and electrons from sources (b) and (c); transverse momentum cut to 
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Figure 4. Left: D° production cross section vs. p t , as it can be measured by 
ALICE in pp (10 9 events) and in central Pb-Pb (10 7 events) collisions; statistical 
(bars) and pt-dependent systematic errors (band) are shown; a normalization error 
of 5% in pp and 11% in Pb-Pb is not shown; pQCD predictions for different sets of 
parameters (m c [GeV], /^/to^c, Mi?/™t,c, PDF set) are also reported for the pp case. 
Right: ratios of simulated ALICE D° data to pQCD curves (parameters: m c [GeV], 
Q/mt,c — HF/mt.c — ^/mt,c, PDF set); the data contain the enhancement due to 
nonlinear gluon evolution while the theory curves do not [9]. 
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Figure 5. Left: ALICE sensitivity for D°-meson Raa with 10 7 central Pb-Pb events 
and 10 9 pp events; statistical (bars) and systematic (bands) errors are shown. Right: 
B production cross section vs. p™ ln reconstructed by ALICE with single muons and 
di-muons in 10 7 central Pb-Pb collisions; only (very small) statistical errors shown. 



reject electrons from charm decays. As an example, with do > 180 fim and p t > 2 GeV, 
the expected statistics of electrons from B decays is 5 x 10 4 for 10 7 central Pb-Pb events, 
with a contamination of about 10%, mainly given by electrons from charm decays [4,24]. 
The sensitivity on the extraction of the bb production cross section and of the B-meson 
p t distribution is currently being investigated. 
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Beauty via (di-)muons in ALICE. B production in Pb-Pb collisions can be measured 
also in the ALICE forward muon spectrometer, 2.5 < rj < 4, analyzing the single-muon 
Pt distribution and the opposite-sign di-muons invariant mass distribution [24]. 

The main backgrounds to the 'beauty muon' signal are 7r ± , K ± and charm decays. 
The cut p t > 1.5 GeV is applied to all reconstructed muons in order to increase the 
signal-to-background ratio. For the opposite-sign di-muons, the residual combinatorial 
background is subtracted using the technique of event-mixing and the resulting 
distribution is subdivided into two samples: the low-mass region, M M + M - < 5 GeV, 
dominated by muons originating from a single b quark decay through b — > c(— > /i + )/U~ 
(BD same ), and the high-mass region, 5 < M^+^ < 20 GeV, dominated by bb — > 
with each muon coming from a different quark in the pair (BBdis). Both samples have 
a background from cc — > /i + yT and a fit is done to extract the charm- and beauty- 
component yields. The single-muon p t distribution has three components with different 
slopes: K and 7r, charm, and beauty decays. Also in this case a fit technique allows 
to extract a p t distribution of muons from B decays. From the //-level cross sections 
a Monte-Carlo-based procedure is used to compute B-level cross sections for the data 
sets (low-mass high- mass and p t -binned single-muon distribution), each 

set covering a specific B-meson p t > p™ m region, as preliminarly shown in Fig. El (right). 
Since only minimal cuts are applied, the reported statistical errors are very small and 
high-p t reach is excellent. Systematic errors are currently under study. 

4. Summary 

We have discussed how heavy quarks, abundantly produced at LHC energies, will allow 
to address several physics issues, in proton-proton, proton-nucleus and nucleus-nucleus 
collisions. In particular, they provide tools to: 

• probe, via parton energy loss and its predicted mass dependence, the high-density 
QCD medium formed in Pb-Pb collisions; 

• probe, in pp collisions, the pQCD calculations parameters space; 

• probe the small-x regime of the PDFs, where saturation/recombination effects are 
expected to be important, even in pp collisions. 

The excellent tracking, vertexing and particle identification performance of ALICE, 
ATLAS and CMS will allow to fully explore this rich phenomenology, as we have shown 
with some specific ALICE studies on D and B meson measurements. 
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